A general exposome factor explains individual differences in functional brain network topography and cognition in youth.

Childhood environments are critical in shaping cognitive neurodevelopment. With the increasing availability of large-scale neuroimaging datasets with deep phenotyping of childhood environments, we can now build upon prior studies that have considered relationships between one or a handful of environmental and neuroimaging features at a time. Here, we characterize the combined effects of hundreds of inter-connected and co-occurring features of a child's environment ("exposome") and investigate associations with each child's unique, multidimensional pattern of functional brain network organization ("functional topography") and cognition. We apply data-driven computational models to measure the exposome and define personalized functional brain networks in pre-registered analyses. Across matched discovery (n=5139, 48.5% female) and replication (n=5137, 47.1% female) samples from the Adolescent Brain Cognitive Development study, the exposome was associated with current (ages 9-10) and future (ages 11-12) cognition. Changes in the exposome were also associated with changes in cognition after accounting for baseline scores. Cross-validated ridge regressions revealed that the exposome is reflected in functional topography and can predict performance across cognitive domains. Importantly, a single measure capturing a child's exposome could more accurately and parsimoniously predict cognition than a wealth of personalized neuroimaging data, highlighting the importance of children's complex, multidimensional environments in cognitive neurodevelopment.

Functional connectivity development along the sensorimotor-association axis enhances the cortical hierarchy.

Human cortical maturation has been posited to be organized along the sensorimotor-association axis, a hierarchical axis of brain organization that spans from unimodal sensorimotor cortices to transmodal association cortices. Here, we investigate the hypothesis that the development of functional connectivity during childhood through adolescence conforms to the cortical hierarchy defined by the sensorimotor-association axis. We tested this pre-registered hypothesis in four large-scale, independent datasets (total n = 3355; ages 5-23 years): the Philadelphia Neurodevelopmental Cohort (n = 1207), Nathan Kline Institute-Rockland Sample (n = 397), Human Connectome Project: Development (n = 625), and Healthy Brain Network (n = 1126). Across datasets, the development of functional connectivity systematically varied along the sensorimotor-association axis. Connectivity in sensorimotor regions increased, whereas connectivity in association cortices declined, refining and reinforcing the cortical hierarchy. These consistent and generalizable results establish that the sensorimotor-association axis of cortical organization encodes the dominant pattern of functional connectivity development.

Generalizable Links Between Borderline Personality Traits and Functional Connectivity.

BACKGROUND: Symptoms of borderline personality disorder (BPD) often manifest during adolescence, but the underlying relationship between these debilitating symptoms and the development of functional brain networks is not well understood. Here, we aimed to investigate how multivariate patterns of functional connectivity are associated with borderline personality traits in large samples of young adults and adolescents. METHODS: We used functional magnetic resonance imaging data from young adults and adolescents from the HCP-YA (Human Connectome Project Young Adult) (n = 870, ages 22-37 years, 457 female) and the HCP-D (Human Connectome Project Development) (n = 223, ages 16-21 years, 121 female). A previously validated BPD proxy score was derived from the NEO Five-Factor Inventory. A ridge regression model with cross-validation and nested hyperparameter tuning was trained and tested in HCP-YA to predict BPD scores in unseen data from regional functional connectivity. The trained model was further tested on data from HCP-D without further tuning. Finally, we tested how the connectivity patterns associated with BPD aligned with age-related changes in connectivity. RESULTS: Multivariate functional connectivity patterns significantly predicted out-of-sample BPD scores in unseen data in young adults (HCP-YA ppermuted = .001) and older adolescents (HCP-D ppermuted = .001). Regional predictive capacity was heterogeneous; the most predictive regions were found in functional systems relevant for emotion regulation and executive function, including the ventral attention network. Finally, regional functional connectivity patterns that predicted BPD scores aligned with those associated with development in youth. CONCLUSIONS: Individual differences in functional connectivity in developmentally sensitive regions are associated with borderline personality traits.

Heritability estimation of subcortical volumes in a multi-ethnic multi-site cohort study.

Heritability of regional subcortical brain volumes (rSBVs) describes the role of genetics in middle and inner brain development. rSBVs are highly heritable in adults but are not characterized well in adolescents. The Adolescent Brain Cognitive Development study (ABCD), taken over 22 US sites, provides data to characterize the heritability of subcortical structures in adolescence. In ABCD, site-specific effects co-occur with genetic effects which can bias heritability estimates. Existing methods adjusting for site effects require additional steps to adjust for site effects and can lead to inconsistent estimation. We propose a random-effect model-based method of moments approach that is a single step estimator and is a theoretically consistent estimator even when sites are imbalanced and performs well under simulations. We compare methods on rSBVs from ABCD. The proposed approach yielded heritability estimates similar to previous results derived from single-site studies. The cerebellum cortex and hippocampus were the most heritable regions (> 50%).

Neural Circuit Markers of Familial Risk for Depression Among Healthy Youths in the Adolescent Brain Cognitive Development Study.

BACKGROUND: Family history of depression is a robust predictor of early-onset depression, which may confer risk through alterations in neural circuits that have been implicated in reward and emotional processing. These alterations may be evident in youths who are at familial risk for depression but who do not currently have depression. However, the identification of robust and replicable findings has been hindered by few studies and small sample sizes. In the current study, we sought to identify functional connectivity (FC) patterns associated with familial risk for depression. METHODS: Participants included healthy (i.e., no lifetime psychiatric diagnoses) youths at high familial risk for depression (HR) (n = 754; at least one parent with a history of depression) and healthy youths at low familial risk for psychiatric problems (LR) (n = 1745; no parental history of psychopathology) who were 9 to 10 years of age and from the Adolescent Brain Cognitive Development (ABCD) Study sample. We conducted whole-brain seed-to-voxel analyses to examine group differences in resting-state FC with the amygdala, caudate, nucleus accumbens, and putamen. We hypothesized that HR youths would exhibit global amygdala hyperconnectivity and striatal hypoconnectivity patterns primarily driven by maternal risk. RESULTS: HR youths exhibited weaker caudate-angular gyrus FC than LR youths (α = 0.04, Cohen's d = 0.17). HR youths with a history of maternal depression specifically exhibited weaker caudate-angular gyrus FC (α = 0.03, Cohen's d = 0.19) as well as weaker caudate-dorsolateral prefrontal cortex FC (α = 0.04, Cohen's d = 0.21) than LR youths. CONCLUSIONS: Weaker striatal connectivity may be related to heightened familial risk for depression, primarily driven by maternal history. Identifying brain-based markers of depression risk in youths can inform approaches to improving early detection, diagnosis, and treatment.

XCP-D: A Robust Pipeline for the post-processing of fMRI data.

Functional neuroimaging is an essential tool for neuroscience research. Pre-processing pipelines produce standardized, minimally pre-processed data to support a range of potential analyses. However, post-processing is not similarly standardized. While several options for post-processing exist, they tend not to support output from disparate pre-processing pipelines, may have limited documentation, and may not follow BIDS best practices. Here we present XCP-D, which presents a solution to these issues. XCP-D is a collaborative effort between PennLINC at the University of Pennsylvania and the DCAN lab at the University at Minnesota. XCP-D uses an open development model on GitHub and incorporates continuous integration testing; it is distributed as a Docker container or Singularity image. XCP-D generates denoised BOLD images and functional derivatives from resting-state data in either NifTI or CIFTI files, following pre-processing with fMRIPrep, HCP, and ABCD-BIDS pipelines. Even prior to its official release, XCP-D has been downloaded >3,000 times from DockerHub. Together, XCP-D facilitates robust, scalable, and reproducible post-processing of fMRI data.

Personalized functional brain network topography is associated with individual differences in youth cognition.

Individual differences in cognition during childhood are associated with important social, physical, and mental health outcomes in adolescence and adulthood. Given that cortical surface arealization during development reflects the brain's functional prioritization, quantifying variation in the topography of functional brain networks across the developing cortex may provide insight regarding individual differences in cognition. We test this idea by defining personalized functional networks (PFNs) that account for interindividual heterogeneity in functional brain network topography in 9-10 year olds from the Adolescent Brain Cognitive Development℠ Study. Across matched discovery (n = 3525) and replication (n = 3447) samples, the total cortical representation of fronto-parietal PFNs positively correlates with general cognition. Cross-validated ridge regressions trained on PFN topography predict cognition in unseen data across domains, with prediction accuracy increasing along the cortex's sensorimotor-association organizational axis. These results establish that functional network topography heterogeneity is associated with individual differences in cognition before the critical transition into adolescence.

Generalizable links between symptoms of borderline personality disorder and functional connectivity.

Background |: Symptoms of borderline personality disorder (BPD) often manifest in adolescence, yet the underlying relationship between these debilitating symptoms and the development of functional brain networks is not well understood. Here we aimed to investigate how multivariate patterns of functional connectivity are associated with symptoms of BPD in a large sample of young adults and adolescents. Methods |: We used high-quality functional Magnetic Resonance Imaging (fMRI) data from young adults from the Human Connectome Project: Young Adults (HCP-YA; N = 870, ages 22-37 years, 457 female) and youth from the Human Connectome Project: Development (HCP-D; N = 223, age range 16-21 years, 121 female). A previously validated BPD proxy score was derived from the NEO Five Factor Inventory (NEO-FFI). A ridge regression model with 10-fold cross-validation and nested hyperparameter tuning was trained and tested in HCP-YA to predict BPD scores in unseen data from regional functional connectivity, while controlling for in-scanner motion, age, and sex. The trained model was further tested on data from HCP-D without further tuning. Finally, we tested how the connectivity patterns associated with BPD aligned with age-related changes in connectivity. Results |: Multivariate functional connectivity patterns significantly predicted out-of-sample BPD proxy scores in unseen data in both young adults (HCP-YA; pperm = 0.001) and older adolescents (HCP-D; pperm = 0.001). Predictive capacity of regions was heterogeneous; the most predictive regions were found in functional systems relevant for emotion regulation and executive function, including the ventral attention network. Finally, regional functional connectivity patterns that predicted BPD proxy scores aligned with those associated with development in youth. Conclusion |: Individual differences in functional connectivity in developmentally-sensitive regions are associated with the symptoms of BPD.

A reproducible and generalizable software workflow for analysis of large-scale neuroimaging data collections using BIDS Apps.

Neuroimaging research faces a crisis of reproducibility. With massive sample sizes and greater data complexity, this problem becomes more acute. Software that operates on imaging data defined using the Brain Imaging Data Structure (BIDS) - BIDS Apps - have provided a substantial advance. However, even using BIDS Apps, a full audit trail of data processing is a necessary prerequisite for fully reproducible research. Obtaining a faithful record of the audit trail is challenging - especially for large datasets. Recently, the FAIRly big framework was introduced as a way to facilitate reproducible processing of large-scale data by leveraging DataLad - a version control system for data management. However, the current implementation of this framework was more of a proof of concept, and could not be immediately reused by other investigators for different use cases. Here we introduce the BIDS App Bootstrap (BABS), a user-friendly and generalizable Python package for reproducible image processing at scale. BABS facilitates the reproducible application of BIDS Apps to large-scale datasets. Leveraging DataLad and the FAIRly big framework, BABS tracks the full audit trail of data processing in a scalable way by automatically preparing all scripts necessary for data processing and version tracking on high performance computing (HPC) systems. Currently, BABS supports jobs submissions and audits on Sun Grid Engine (SGE) and Slurm HPCs with a parsimonious set of programs. To demonstrate its scalability, we applied BABS to data from the Healthy Brain Network (HBN; n=2,565). Taken together, BABS allows reproducible and scalable image processing and is broadly extensible via an open-source development model.

Polyneuro risk scores capture widely distributed connectivity patterns of cognition.

Resting-state functional connectivity (RSFC) is a powerful tool for characterizing brain changes, but it has yet to reliably predict higher-order cognition. This may be attributed to small effect sizes of such brain-behavior relationships, which can lead to underpowered, variable results when utilizing typical sample sizes (N∼25). Inspired by techniques in genomics, we implement the polyneuro risk score (PNRS) framework - the application of multivariate techniques to RSFC data and validation in an independent sample. Utilizing the Adolescent Brain Cognitive Development® cohort split into two datasets, we explore the framework's ability to reliably capture brain-behavior relationships across 3 cognitive scores - general ability, executive function, learning & memory. The weight and significance of each connection is assessed in the first dataset, and a PNRS is calculated for each participant in the second. Results support the PNRS framework as a suitable methodology to inspect the distribution of connections contributing towards behavior, with explained variance ranging from 1.0 % to 21.4 %. For the outcomes assessed, the framework reveals globally distributed, rather than localized, patterns of predictive connections. Larger samples are likely necessary to systematically identify the specific connections contributing towards complex outcomes. The PNRS framework could be applied translationally to identify neurologically distinct subtypes of neurodevelopmental disorders.

Development of top-down cortical propagations in youth.

Hierarchical processing requires activity propagating between higher- and lower-order cortical areas. However, functional neuroimaging studies have chiefly quantified fluctuations within regions over time rather than propagations occurring over space. Here, we leverage advances in neuroimaging and computer vision to track cortical activity propagations in a large sample of youth (n = 388). We delineate cortical propagations that systematically ascend and descend a cortical hierarchy in all individuals in our developmental cohort, as well as in an independent dataset of densely sampled adults. Further, we demonstrate that top-down, descending hierarchical propagations become more prevalent with greater demands for cognitive control as well as with development in youth. These findings emphasize that hierarchical processing is reflected in the directionality of propagating cortical activity and suggest top-down propagations as a potential mechanism of neurocognitive maturation in youth.

Curation of BIDS (CuBIDS): A workflow and software package for streamlining reproducible curation of large BIDS datasets.

The Brain Imaging Data Structure (BIDS) is a specification accompanied by a software ecosystem that was designed to create reproducible and automated workflows for processing neuroimaging data. BIDS Apps flexibly build workflows based on the metadata detected in a dataset. However, even BIDS valid metadata can include incorrect values or omissions that result in inconsistent processing across sessions. Additionally, in large-scale, heterogeneous neuroimaging datasets, hidden variability in metadata is difficult to detect and classify. To address these challenges, we created a Python-based software package titled "Curation of BIDS" (CuBIDS), which provides an intuitive workflow that helps users validate and manage the curation of their neuroimaging datasets. CuBIDS includes a robust implementation of BIDS validation that scales to large samples and incorporates DataLad--a version control software package for data--as an optional dependency to ensure reproducibility and provenance tracking throughout the entire curation process. CuBIDS provides tools to help users perform quality control on their images' metadata and identify unique combinations of imaging parameters. Users can then execute BIDS Apps on a subset of participants that represent the full range of acquisition parameters that are present, accelerating pipeline testing on large datasets.

Learning Strategies for Contrast-agnostic Segmentation via SynthSeg for Infant MRI data.

Longitudinal studies of infants' brains are essential for research and clinical detection of neurodevelopmental disorders. However, for infant brain MRI scans, effective deep learning-based segmentation frameworks exist only within small age intervals due to the large image intensity and contrast changes that take place in the early postnatal stages of development. However, using different segmentation frameworks or models at different age intervals within the same longitudinal data set would cause segmentation inconsistencies and age-specific biases. Thus, an age-agnostic segmentation model for infants' brains is needed. In this paper, we present "Infant-SynthSeg", an extension of the contrast-agnostic SynthSeg segmentation framework applicable to MRI data of infants at ages within the first year of life. Our work mainly focuses on extending learning strategies related to synthetic data generation and augmentation, with the aim of creating a method that employs training data capturing features unique to infants' brains during this early-stage development. Comparison across different learning strategy settings, as well as a more-traditional contrast-aware deep learning model (nnU-net) are presented. Our experiments show that our trained Infant-SynthSeg models show consistently high segmentation performance on MRI scans of infant brains throughout the first year of life. Furthermore, as the model is trained on ground truth labels at different ages, even labels that are not present at certain ages (such as cerebellar white matter at 1 month) can be appropriately segmented via Infant-SynthSeg across the whole age range. Finally, while Infant-SynthSeg shows consistent segmentation performance across the first year of life, it is outperformed by age-specific deep learning models trained for a specific narrow age range.